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> ■ ABSTRACT 

OO I Successfully reproducing the galaxy luminosity function and the bimodality in the 

i galaxy distribution requires a mechanism that can truncate star formation in massive 

\l ' haloes. Current models of galaxy formation consider two such truncation mechanisms: 

strangulation, which acts on satellite galaxies, and AGN feedback, which predomi- 
nantly affects central galaxies. The efficiencies of these processes set the blue fraction 
of galaxies, fhiue{L, M), as function of galaxy luminosity, L, and halo mass, M. In this 
paper we use a galaxy group catalogue extracted from the Sloan Digital Sky Survey 
• (SDSS) to determine fhinc{L, M). To demonstrate the potential power of these data 

as a benchmark for galaxy formation models, we compare the results to the semi- 
Q . analytical model for galaxy formation of Croton et al. (2006). Although this model 

' accurately fits the global statistics of the galaxy population, as well as the shape of 

I the conditional luminosity function, there are significant discrepancies when the blue 

■ fraction of galaxies as a function of mass and luminosity is compared between the ob- 

servations and the model. In particular, the model predicts (i) too many faint satellite 
galaxies in massive haloes, (ii) a blue fraction of satellites that is much too low, and 
(iii) a blue fraction of centrals that is too high and with an inverted luminosity de- 
I pendence. In the same order, we argue that these discrepancies owe to (i) the neglect 

of tidal stripping in the semi-analytical model, (ii) the oversimplified treatment of 
strangulation, and (iii) improper modeling of dust extinction and/or AGN feedback. 
The data presented here will prove useful to test and calibrate future models of galaxy 
formation and in particular to discriminate between various models for AGN feedback 
and other star formation truncation mechanisms. 
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1 INTRODUCTION 

One of the most challenging outstanding problems in galaxy 
formation is to explain the detailed shape of the galaxy 
luminosity function (hereafter LF). In particular, the rel- 
atively shallow faint-end slope and the exponential cut-off 
at the bright end of the LF have proven difficult to explain 
(e.g.. White & Frenk 1991; Benson et al. 2003). In the tra- 
ditional scenario for galaxy formation, it is envisioned that 
lower cooling efficiencies in massive galaxies would explain 
the exponential tail of the LF (Rees & Ostriker 1977; Silk 
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1977; White & Rees 1978), while supernova feedback is typ- 
ically invoked to reduce the star formation efficiency in low 
mass haloes (Larson 1974; White & Rees 1978; Dekel & Silk 
1986). Although the latter can indeed be tuned to reproduce 
the faint-end slope of the galaxy LF, it worsens the prob- 
lems at the bright end. As nicely demonstrated by Benson 
et al. (2003), supernova feedback causes a drastic increase 
of the amount of diffuse hot gas that remains in larger ha- 
los. This gas is able to cool onto the central galaxies in these 
haloes, producing too many bright galaxies. In addition, this 
causes the bright model galaxies to have relatively young 
stellar populations, in disagreement with observations (e.g., 
Kauffmann & Chariot 1998; Heavens et al. 2004; Thomas 
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et al. 2005). What is needed is a mechanism that can trun- 
cate the star formation in these massive, central galaxies at 
relatively late times. 

Star formation truncation is also the main mechanism 
that is thought to underlie the bimodality of galaxy prop- 
erties. The local population of galaxies consists roughly of 
two types: red galaxies, which reveal an early type morphol- 
ogy and which have very little or no ongoing star formation, 
and blue galaxies with active star formation and a late-type 
morphology (e.g., Stratcva ct al. 2001; Blanton et al. 2003a; 
Kauffmann et al. 2003, 2004; Baldry et al. 2004; Brinchmann 
et al. 2004; Balogh et al. 2004a,b). Although a non-negligible 
fraction of red galaxies arc clearly edge-on disk galaxies that 
owe their red colour to an enhanced extinction, the most 
pronounced distinction between 'red sequence' and 'blue se- 
quence' galaxies is their current star formation rate. Since 
only relatively small amounts of ongoing star formation tend 
to make a galaxy appear 'blue', the colour bimodality basi- 
cally reflects star formation truncation: red sequence galax- 
ies have their star formation truncated, while blue sequence 
galaxies are still forming stars today. 

Semi-analytical models for galaxy formation consider 
a number of mechanisms that can prevent, delay or trun- 
cate star formation. In low mass haloes one typically invokes 
reionization and supernova feedback in order to suppress or 
truncate star formation. In haloes with M >, W^^h^^ M©, 
on which we will focus in this paper, two additional trunca- 
tion effects are considered. The first one, called strangula- 
tion, only aifects satellite galaxies. As soon as a dark mat- 
ter halo is accreted by a larger halo, its central galaxy be- 
comes a satellite galaxy. It is often assumed that this accre- 
tion process causes the satellite galaxy to be stripped of its 
hot gas reservoir. Consequently, after a delay time in which 
the galaxy consumes (part of) its cold gas, star formation 
is truncated, and the satellite galaxy becomes red (Larson, 
Tinslcy & Caldwell 1980; Balogh, Navarro & Morris 2000). 
Ram-pressure stripping (Gunn & Gott 1972) may shorten 
the time-delay, by also stripping the satellite of its cold gas 
reservoir, but since the time scale for strangulation is already 
relatively short, the addition of ram- pressure stripping does 
not have a large impact. Virtually all semi-analytical mod- 
els of galaxy formation, starting with Kauffmann, White 
& Guiderdoni (1993), have talcen this strangulation mech- 
anism into account. In fact, it is the main mechanism that 
causes red-sequence model galaxies to preferentially reside 
in overdcnsc regions such as groups and clusters of galax- 
ies, in good agreement with observations (e.g. Oemler 1974; 
Dressier 1980; Hogg et al. 2004; Balogh et al. 2004b; Wein- 
mann et al. 2006). 

The second star formation truncation mechanism that 
operates in massive haloes is feedback from active galactic 
nuclei (AGN), which predominantly affects central galaxies. 
Although the potential importance of AGN feedback has 
long been recognized (e.g., Tabor & Binney 1993; Ciotti & 
Ostriker 1997), it only recently has been given serious con- 
sideration in galaxy formation models. This has largely been 
motivated by X-ray observations which reveal that AGN can 
indeed impact the hot IGM of galaxy clusters (e.g., Fabian 
ct al. 2003; McNamara ct al. 2005). Numerous recent stud- 
ies have demonstrated that the inclusion of AGN feedback 
in galaxy formation models can help to explain the bright 
end exponential cut-off of the galaxy luminosity function 



(e.g. Granato et al. 2004; Croton et al. 2006; Cattaneo et 
al. 2006a; Sijacki & Springel 2005; Bower et al. 2005) and 
the fact that the most massive galaxies contain the oldest 
stars (Croton et al. 2006; Scannapieco, Silk & Bouwens 2005; 
Bower et al. 2005). 

Despite these successes, we are still far from a proper 
understanding of how AGN feedback may establish an equi- 
librium state where it can efficiently suppress star formation 
in the centers of massive haloes. In the studies mentioned 
above, AGN feedback is typically modeled using oversimpli- 
fied, heuristic scaling relations, often based on very different 
views of how AGN feedback might operate. In particular, it 
is still unclear which mode of AGN activity is most impor- 
tant for the star formation truncation discussed above; the 
merger induced 'quasar mode' which leads to an initial star- 
burst followed by a quenching of star formation, or the 'radio 
mode', which is caused by continual and quiescent accretion 
of hot gas onto the central supcrmassive black hole. Hop- 
kins et al. (2006) have suggested that merger-induced AGN 
activity (the 'quasar mode') is responsible for the transition 
from blue, star forming to red, passive galaxies. Springel, Di 
Matteo & Hernquist (2005a), Menci et al. (2005) and Kang, 
Jing & Silk (2005) have shown that this 'quasar mode' feed- 
back can indeed terminate star formation and expel the gas 
from the center of the galaxy, once the supermassive black 
holes become sufficiently massive. On the other hand, Cro- 
ton et al. (2006), Bower et al. (2005), Cattaneo et al. (2006a), 
Nusser ct al. (2006) and Sijacki & Springel (2006) have ar- 
gued that the 'radio mode' of AGN activity (or AGN feed- 
back operating in quasi-hydrostatically cooling haloes) is 
the main mechanism to truncate star formation in massive 
galaxies. How exactly this radio mode feedback operates, 
however, is still unclear, as is evident from the fact that the 
aforementioned studies all use very different formulations. 

All these different AGN feedback models mainly differ 
in the way in which the feedback efficiency scales with halo 
mass and with galaxy properties such as black hole mass and 
gas mass fraction. Since AGN feedback causes star formation 
truncation, one way to discriminate between these various 
models is therefore to investigate the relative fractions of 
blue and red galaxies as function of halo mass and galaxy 
properties. Such a study will also help to improve our un- 
derstanding of strangulation, the stax formation truncation 
mechanism for satellites. Although strangulation is typically 
modeled as being independent of halo mass, one might ar- 
gue that the ability of a host halo to strip a subhalo of its 
hot gas reservoir depends on the presence and density of the 
hot corona of the host halo, which in turn may well be mass 
dependent. Again, knowledge of the fractions of blue and 
red (satellite) galaxies as function of halo mass should allow 
us to discriminate between these different possibilities. 

Nowadays, with large galaxy surveys, such as the Two 
Degree Field Galaxy Redshift Survey (Colless et al. 2001) 
and the Sloan Digital Sky Survey (York et al. 2000), the 
number of galaxies is sufficiently large that, in principle, 
one could accurately measure the red and blue fractions as 
a function of various variables. In this paper, we use our 
SDSS group catalogue, presented in Weinmann et al. (2006, 
hereafter paper I), to compute the fractions of blue and 
red galaxies as function of both halo mass and galaxy lu- 
minosity. To emphasize the potential constraining power of 
these data we compare our results to the semi-analytical 
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galaxy formation model of Croton ct al. (2006), which in- 
cludes both strangulation and 'radio-mode' AGN feedback. 
We show that although this model accurately fits the galaxy 
luminosity function, the color-bimodality, and many other 
global statistics of the galaxy population, it fails dramati- 
cally when it comes down to the blue fraction of galax;ies 
as a function of halo mass and luminosity. We argue that 
this has its origin in the way that strangulation and AGN 
feedback have been incorporated and we briefly discuss pos- 
sible modifications. The aim of this paper, however, is not 
to present a new, improved model for star formation trun- 
cation, but merely to present observational constraints that 
will hopefully prove useful in discriminating between the 
various models. 

This paper is organized as follows. Section 2 describes 
our SDSS group catalogue, which we compare to a similar 
group catalogue extracted from the semi-analytical model of 
Croton et al. (2006) described in Section 3. The actual com- 
parison is presented in Section 4, while Section 5 discusses 
the possible implications for galaxy formation models. We 
summarize our results in Section 6. 



2 THE SDSS GROUP CATALOGUE 

2.1 Data 

The Sloan Digital Sky Survey (SDSS; York ct al. 2000) 
is a joint, five passband («, g, r, i and z) imaging and 
medium-resolution [R ^ 1800) spectroscopic survey. In this 
paper, we focus on the subset of galaxies that are in the 
New York University Value-Added Galaxy Catalogue (NYU- 
VAGC) based on the SDSS Data Release 2 (Blanton et 
al. 2005). This NYU-VAGC is based on an independent, 
significantly improved data reduction. From this catalogue 
we select all galaxies with an extinction corrected apparent 
magnitude brighter than r = 17.77, with redshifts in the 
range 0.01 < z < 0.20, and with a rcdshift completeness 
C > 0.7. This leaves a grand total of 184,425 galaxies. In 
what follows, we use Mr and '^'^ Mr to indicate the absolute 
magnitude in the r-band, k-corrected to 2: = and z = 0.1, 
respectively. All k-corrections are based on the model de- 
scribed in Blanton et al. (2003b). 

We split the galaxies into 'red' and 'blue' subsamples 
using a magnitude dependent cut, which roughly follows the 
observed bimodality scale in the colour-magnitude relation: 

°'\9 - Ocut = 0.7 - 0.032 ["-^M, + 16.5] (1) 

(cf. paper I). In what follows, we refer to galaxies that are 
redder and bluer than — r)cut as 'red' and 'blue' galax- 
ies, respectively. 

2.2 The group-finding algorithm 

Our working definition of a galaxy group is the ensemble of 
galaxies that reside in the same dark matter parent halo; 
galaxies that reside in subhaloes are considered to be group 
members that belong to the parent halo in which the sub- 
halo is located. We have used the halo-based group finder 
developed by Yang et al. (2005a, hereafter YMBJ) in order 
to group the galaxies in the above mentioned galaxy cat- 
alogue. This particular group finder has been optimized to 



group galaxies according to their common dark matter halo, 
and has been thoroughly tested with mock galaxy redshift 
surveys. In brief, the method works as follows. First, poten- 
tial group centers are identified using a Friends-Of-Friends 
(POP) algorithm or an isolation criterion. Next, the total 
group luminosity is estimated which is converted into an es- 
timate for the group mass using an assumed mass-to-light 
ratio. From this mass estimate, the radius and velocity dis- 
persion of the corresponding dark matter halo are estimated 
using the virial equations, which in turn are used to select 
group members in redshift space. This method is iterated 
until group memberships converge. A more detailed descrip- 
tion is given in Appendix A of paper I. 

In YMB.J the performance of this group finder has been 
tested in terms of completeness of true members and con- 
tamination by interlopers, using detailed mock galaxy red- 
shift surveys. The average completeness of individual groups 
was found to be ~ 90 percent, with only ~ 20 percent inter- 
lopers. Furthermore, the resulting group catalogue is insensi- 
tive to the initial assumption regarding the mass-to-light ra^- 
tios, and the group finder is more successful than the conven- 
tional FOF method (e.g., Huchra &: Geller 1982; Ramella, 
GcUcr & Huchra 1989; Merchan & Zandivarcz 2002; Eke ct 
al. 2004; Berlind et al. 2006) in associating galaxies accord- 
ing to their common dark matter haloes. 

2.3 Estimating group masses 

Following YMBJ we use the group luminosity to assign 
masses to our groups. The motivation behind this is that 
one naturally expects the group luminosity to be strongly 
correlated with halo mass (albeit with a certain amount of 
scatter). Since the group luminosity is dominated by the 
brightest members, which are exactly the ones that can be 
observed in a flux limited survey like the SDSS, the determi- 
nation of the (total) group luminosity is more robust than 
that of the group's velocity dispersion, especially when the 
number of group members is small (see Appendix B in pa- 
per I). 

Clearly, because of the flux limit of the SDSS, two iden- 
tical groups observed at diff'ercnt redshifts will have a differ- 
ent Lgroup, defined as the summed luminosity of all its identi- 
fied members. To circumvent this bias we first need to bring 
the group luminosities to a common scale. A nearby group 
selected in an apparent magnitude limited survey should 
contain all of its members down to a faint luminosity. We can 
therefore use these nearby groups to determine the relation 
between the group luminosity obtained using only galax- 
ies above a bright luminosity limit and that obtained using 
galaxies above a fainter luminosity limit. Assuming that this 
relation is redshift-independent, one can correct the lumi- 
nosity of a high-0 group, where only the brightest members 
are observed, to an empirically normalized luminosity scale. 

As common luminosity scale we use Lia.r, , defined as 
the luminosity of all group members brighter than "'^Mr = 

— 19.5 + 5 log /i. To calibrate the relation between Lgroup and 
Lia.5 we first select all groups with z < 0.09, which cor- 
responds to the redshift for which a galaxy with "'^Mr = 

— 19.5 + 5 log /i has an apparent magnitude that is equal to 
the magnitude limit of the survey. For groups with z > 0.09 
we use this 'local' calibration between Lgroup and Lig.s to es- 
timate the latter. Detailed tests have shown that the result- 
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ing group luminosities are significantly more reliable than 
those in which the correction for missing members is based 
on the assumption of a universal luminosity function (see 
YMBJ for details). 

The final step is to obtain an estimate of the group 
(halo) mass from Lig.s. This is done using the assumption 
that there is a one-to-one relation between L19.5 and halo 
mass. For each group we determine the number density of 
all groups brighter (in terms of Lig.s) than the group in 
consideration. Using the halo mass function corresponding 
to a ACDM concordance cosmology with flm = 0.3, SIa = 
0.7, h = Ho/(100 kms"^Mpc"^) = 0.7 and erg = 0.9 we 
then find the mass for which the more massive haloes have 
the same number density. Although the masses thus derived 
depend on cosmology, it is straightforward to convert the 
masses derived here to any other cosmology. 

Finally we note that not all groups can have a halo mass 
assigned to them. First of all, the mass estimator described 
above does not work for groups in which all members arc 
fainter than '^'^Mr = —19.5 + 5logh. Secondly, the combi- 
nation of 1/19.5 and redshift may be such that we know that 
the halo catalogue is incomplete, which means that there 
is a significant rmmber of groups at this redshift with the 
same L19.5 but for which the individual galaxies are too 
faint to be detected. Since our mass assignment is based on 
the assumption of completeness, any group beyond the com- 
pleteness redshift corresponding to its L19.5 is not assigned 
a halo mass (see Yang et al. 2005b for details). 

2.4 The SDSS group catalogue 

Applying our group finder to the sample of SDSS galaxies 
described in Section 2.1 yields a group catalogue of 53,229 

systems with an estimated mass. Those groups contain a 
total of 92,315 galaxies. The majority of the groups (37,216 
systems) contain only a single member, while there are 9220 
binary systems, 3073 triplet systems, and 3720 systems with 
four members or more (see paper I for details)^. In what 
follows we refer to the brightest galaxy in each group as the 
'central' galaxy, while all others are termed 'satellites'. 



3 THE SEMI-ANALYTIC MODEL 

The semi-analytic model (hereafter SAM) to which we com- 
pare the SDSS group catalogue discussed above is based 
on an output of the Millennium Run A'^-body simulation 
(Springel et al. 2005b) and is described in detail in Cro- 
ton et al. (2006; hereafter C06). The simulation is based 
on the cosmological parameters Qm = 0.25, f^A = 0.75, 
Qh = 0.045, h = 0.73 and erg = 0.9 and has a volume of 
0.125/i~^Gpc^. Dark matter haloes are identified with a FOF 
group finder, and subsequently populated with galaxies fol- 
lowing the semi-analytical model described in C06. 

One of the relative novelties of this SAM is the inclusion 
of 'radio mode' feedback from AGN that lie at the center of 
a halo with a static corona of hot gas. As shown in C06, 
this feedback mode suppresses the cooling flow in massive 

^ This SDSS group catalogue is publicly available at 
http: / /www.astro.umass.edu/^xhyang/Group.html 



haloes at relatively late times, which in turn yields luminosi- 
ties, colours and stellar ages for massive galaxies in better 
agreement with observations. In particular, the inclusion of 
the radio-mode AGN feedback can explain the exponential 
cut-off at the bright end of the galaxy luminosity function, 
and the fact that the most massive galax;ies are red and 
consist of old stellar populations (see De Lucia et al. 2006 
for details). The model also predicts stax formation histo- 
ries, cold gas mass fractions and metallicities that are all in 
good agreement with observations. The model even predicts 
a TuUy-Fisher zero-point that matches the data, as long as 
the rotation velocity of a disk galaxy is equal to the maxi- 
mum circular velocity of the dark matter halo (however, see 
Dutton et al. 2006 for the unorthodox implications of this 
assumption) . 

Since the halo masses M assigned to our SDSS galaxy 
groups are obtained by matching the abundances to the halo 
mass function, these are to be interpreted as the masses 
inside a radius with an overdensity of 180. As shown by 
Jenkins et al. (2001), for this definition of halo mass the 
analytical halo mass function of Sheth, Mo & Tormen (2001) 
used here is in good agreement with the halo mass function 
obtained from numerical simulations. The halo masses in the 
C06 catalogue, however, are defined as the masses inside a 
radius with a mean density that is 200 times the critical 
density, which we denote by M200 • In order to convert M200 
to M, we assume that dark matter haloes follow a NFW 
density distribution (Navarro, Frenk & White 1997). Using 
the relation between halo mass and halo concentration of 
Eke, Navarro & Steinmetz (2001), we find that the relation 
between M and M200 is well fit by 

= 0.745 - 0.0006 [log(M2oo) - 7.0]^-*^ . (2) 

In what follows, wo only consider the galaxies with 
Mr < —16.72, which refiects the magnitude completeness 
limit of the SAM, yielding a total of approximately 9 mil- 
lion model galaxies. The main information used in this paper 
are the r band and g band magnitudes of these galaxies, and 
the mass M of the halo in which they reside. For comparison 
with the SDSS, we also compute the g and r-band magni- 
tudes k-corrected to z = 0.1, using 

°-^g = g + 0.3113 + 0.4620 {g-r- 0.6102) (3) 
and 

°V = g- 0.4075 - 0.8577 {g-r- 0.6102) (4) 

(Blanton et al. 2003b)^. 

As we will see below, despite the AGN feedback the 
SAM still contains a significant number of very bright and 
blue galaxies that are not present in the SDSS. As men- 
tioned in C06, these are mainly ULIRG-type starbursts for 
which the dust treatment of the model is inadequate; with 
a more proper dust model, these galaxies would be much 
more extincted, making them both fainter and redder. In 
order to suppress the impact of these galaxies on our SAM- 
SDSS comparison we remove all galaxies with (<? — r) < 
from both the SAM and the SDSS. In the case of the SAM 

^ These filter transformations are taken from a manuscript in 
preparation by Michael Blanton and Sam Roweis, available at 
http: / / cosmo.nyu.edu/blanton/kcorrect 
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Figure 1. A comparison of global statistics of the SAM-RS (right-hand panels) with the SDSS (left-hand panels). Panels in the first and 
second row show histograms of absolute "-^r-band magnitude and "'^(g — r) colour for both redshift surveys. The contributions from blue 
and red galaxies are indicated by dashed and dotted lines, respectively. Note the good agreement between SDSS and SAM-RS. Panels 
in the third row show histograms of morphological parameters. In the case of the SDSS, we plot the distribution of the concentration 
parameter c, defined as the ratio of the radii containing 90 and 50 percent of the petrosian fiux. In the case of the SAM-RS, we plot the 
distribution of the bulge-to-total mass ratio B/T instead. Although these can not be compared directly, in general a more concentrated 
galaxy (higher c) will also have a larger B/T. Again, the contributions from red and blue galaxies are indicated. Finally, the fourth row 
of panels shows the colour-magnitude relations. The solid line indicates the bimodality scale given by equation (1), which wc use to split 
the population of galaxies in red and blue subpopulations. 
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Figure 2. The average mass-to-light ratio, (M/L19.5) as a func- 
tion of halo (group) mass. Here L19.5 is the total luminosity of 
all galaxies in a halo with °-'^Mr - 51og/i < -19.5. Filled circles 
with errorbars (indicating the 68 % confidence level) show the 
results obtained from the SAM directly, using the original halo 
masses and halo membership. The solid line shows the (M/Lig.s) 
as obtained from the SAM group catalogue, and is thus based on 
the assigned group masses and the assigned group memberships. 
The agreement with the true (M/L19.5} is excellent, indicating 
that our group finder allows an accurate recovery of the average 
relation between mass and light. 



this only affects 0.5 percent of all galaxies, while in the case 
of the SDSS, this fraction is completely negligible. 

We split the SAM model galaxies into 'red' and 'blue' 
subsamples using the same magnitude dependent colour-cut 
as for the SDSS, given by equation (1). In addition, we dis- 
criminate between 'central' galaxies, defined as the galaxy in 
a halo that is closest to the halo center, and 'satellite' galax- 
ies. This differs from the definition used for the group cata- 
logues, where the central galaxy is defined as the brightest 
group member. However, 97.5 percent of all central galax- 
ies in the SAM are also the brightest galaxy in their halo, 
virtually independent of halo mass. We have verified that 
defining central galaxies in the SAM as the brightest halo 
members instead does not have a significant impact on any 
of our results. 



3.1 Constructing a SAM redshift survey 

In order to be able to compare the SAM to the SDSS results 
we need to mimic the construction of a galaxy redshift sur- 
vey. We do so using the following steps. First, we construct 
a large virtual universe by replicating the periodic simula- 
tion box in a stack of 2 x 2 x 2 boxes. This is required in 
order to be able to probe out to sufficiently high redshifts. 
Next we compute the redshift and apparent magnitude of 
each galaxy as seen by a virtual observer located in a cor- 
ner of this virtual universe (who can thus see 7r/2 steradian 
of 'sky'). We mimic the selection criteria of the SDSS dis- 



cussed in Section 2.1 by only selecting those galaxies with 
0.01 < z < 0.2 and with r < 17.77. This leaves us with a 
grand total of 428, 013 model galaxies. In what follows, we 
refer to this sample as the "SAM redshift survey" (hereafter 
SAM-RS). 

Fig. 1 compares a number of statistics of these galaxies 
with those from the SDSS. The upper panels plot the distri- 
bution of absolute magnitudes, "'^Mr, and their contribution 
due to red (dotted curves) and blue (dashed curves) galaxies. 
The agreement is very satisfactory, consistent with the fact 
that the SAM matches the observed luminosity functions 
(see C06). The second row of panels indicate the color distri- 
butions. Once again, the agreement is reasonable, although 
the bimodality in the SAM-RS seems somewhat more pro- 
nounced than in the SDSS, with a somewhat narrower 'red 
peak'. Nevertheless, with blue fractions of 48 and 46 percent 
in the SDSS and SAM-RS, respectively, the overall agree- 
ment is very satisfactory. 

In the third row of panels, we 'compare' two different 
morphology indicators. For the SDSS galaxies we plot the 
distribution of the concentration c, defined as the ratio be- 
tween the radii that contain 90 and 50 percent of the Pet- 
rosian fiux. For the SAM model galaxies, we plot the distri- 
bution of the bulge-to-total (B/T) stellar mass ratio instead. 
Typically, a galaxy with a large B/T will also have a high 
concentration parameter. For both c and B/T there is a very 
significant overlap of red and blue galaxies. Therefore, our 
split in red and blue galaxies does not necessarily correspond 
to a morphological split in early and late-type galaxies, re- 
spectively, even though both are clearly correlated (see Fig. 1 
in paper I). 

Finally, the lower panels of Fig. 1 show scatter plots 
of the color-magnitude relations. The solid line corresponds 
to the bimodality scale given by equation (1). Again, there 
is reasonable overall agreement between the SAM-RS and 
the SDSS, although there are more galaxies with very blue 
colours in the SAM-RS, especially at the bright end. Recall 
that galaxies with {g ~ r) < have already been excluded 
from these plots. As mentioned before, these bright, blue 
galaxies would appear significantly fainter and redder with 
proper dust modeling. Another apparent discrepancy con- 
cerns the red sequence, which at the bright end appears 
significantly tighter for the SAM than for the SDSS, which 
is also apparent from the histograms in the second row of 
panels. 

To summarize, despite some small discrepancies, the 
global, statistical properties of the SAM model galaxies are 
in good agreement with the SDSS. However, this does not 
mean that the SAM also predicts the correct statistics as a 
function of halo mass. This is clearly a much tighter con- 
straint for the model, and, as we argued in Section 1, may 
provide useful insights regarding the halo-mass dependence 
of various physical processes. 

3.2 Construction of the SAM group catalogue 

The main goal of this paper is to compare the fractions of 
red and blue galaxies as a function of both halo mass and 
luminosity in the SAM with those obtained from our SDSS 
group catalogue. In principle we could do so by comparing 
our SDSS group results directly with the SAM. However, 
our group finder, and in particular the algorithm used to as- 
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Figure 3. The conditional luminosity function obtained from the SDSS group catalogue (upper panel) from the SAM directly, using 
true halo masses and true halo members (middle panel) and from the SAM group catalogue, using assigned group masses and assigned 
memberships (lower panels). Results are shown for four different mass bins (masses in Mq), as indicated at the top of each column. 
Luminosities are in the '^•^r-band. The dotted (blue) line marks the contribution from the satellite galaxies, the dashed (red) line the 
contribution from the central galaxies. 



sign the masses, is not perfect. Hence, it is unclear whether 
any discrepancy between SDSS group catalogue and SAM 
reflects an artefact of the group finder, or whether there is 
a true difference in the halo occupation statistics of SDSS 
and SAM. To circumvent this problem, we use the SAM- 
RS described above to construct a "SAM group catalogue" 
(hereafter SAM-GC), using exactly the same group finder 
and mass assignment algorithm as for the SDSS. By com- 
paring this SAM-GC with the SDSS group catalogue we sig- 
nificantly reduce the impact of possible inaccuracies related 
to the group finder, making the comparison more fair. 

Application of our group finder to the SAM-RS de- 
scribed above yields 98,130 groups with an assigned mass, 
which host a total of 206,076 galaxies. This amounts to an 
average of 2.10 members per group, which is significantly 
higher than for the SDSS, where the groups with an assigned 
mass have on average 1.73 members. However, the SDSS is 
not complete; in fact, the average completeness of our SDSS 
sample is 0.88, which largely explains the difference in the 
mean number of members per group. Another reason for 
this discrepancy is that, as we will see, massive haloes in 



the SAM contain too many faint satellites compared to the 
SDSS. 

The fraction of blue galaxies in the SAM-GC is only 29 
percent. Comparing this to the fraction of 46 percent of blue 
galaxies in the SAM-RS indicates that a red galaxy is much 
more likely to be associated with a group than a blue galaxy. 
At first sight this seems a logical consequence of the fact that 
(i) our group catalogue is limited to relatively massive haloes 
with >, 5 X 10^^^"^ M0 and (ii) low mass haloes are more 
likely to host blue galaxies, see e.g. paper 1. However, the 
application of the group finder to the SDSS only reduces the 
fraction of blue galaxies from 48 percent to 41 percent. This 
reduction is much less severe than for the SAM. This is the 
first indication that the SAM and SDSS do not agree well 
when it comes down to details regarding the distribution of 
red and blue galaxies (see Section 4.2 below). 

Since the SAM contains the full halo occupation in- 
formation, we can check whether our group finder has as- 
signed the correct galaxies to the same group, and whether 
the assigned mass is in agreement with the true halo mass. 
We have performed a large number of tests to investigate 
how well the group finder allows us to recover the true re- 
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lations between galaxies and their dark matter haloes. Sev- 
eral of these tests have been described in detail in Yang et 
al. (2005a,b), and show that the average occupation statis- 
tic of dark matter haloes are accurately recovered. However, 
higher-order moments of the occupation statistics, such as 
the scatter around a mean relation, are typically severely 
underestimated by the group catalogue, due to the fact that 
we assume a one-to-one relation between halo mass and halo 
luminosity (with zero scatter) when assigning masses to our 
groups. As an illustration of the accuracy of our group finder. 
Fig. 2 plots the average mciss-to-light ratio of the dark mat- 
ter haloes in the SAM (symbols with errorbars). Here M is 
the halo mass defined according to equation (2), and L is the 
total luminosity in the "'^r band of all galaxies in that halo 
with "-^Mr - 5 log /t < -19.5. The solid line in Fig. 2 shows 
the average mass-to- light ratios obtained from the SAM-GC, 
which agree extremely well with the true {M/L)m- This 
demonstrates that our group finder can accurately recover 
the average relation between mass and light. 

In what follows, whenever we present any result ob- 
tained from the SAM-GC, we will also present the same 
results extracted directly from the SAM (using the true 
halo masses and the true halo memberships). A comparison 
among these results safeguards against potential problems 
with the group finder. 



4 GALAXY ECOLOGY 

4.1 Conditional Luminosity Functions 

We start our SAM-SDSS comparison by focusing on the con- 
ditional luminosity function (CLF), "1>(L|M), which specifies 
the average number of galaxies of luminosity L that reside 
in a halo of mass M (Yang, Mo & van den Bosch 2003; van 
den Bosch, Yang & Mo 2003). The upper panels of Fig. 3 
show the CLFs obtained from the SDSS group catalogue. 
Results are shown for four difi^erent mass bins, as indicated 
at the top of each column. Note that these masses are the 
assigned group masses. The dotted (blue) and dashed (red) 
lines indicate the contributions from the satellite and cen- 
tral galaxies, respectively. The distribution of central gala^x- 
ies is well approximated by a log-normal distribution, con- 
sistent with previous findings (Zheng et al. 2005; Yang et 
al. 2005b). The panels in the middle row of Fig. 3 show the 
CLFs obtained directly from the SAM; here the masses are 
the true halo masses, and the true halo members are used to 
construct the CLFs. The overall agreement with the CLFs 
extracted from the SDSS group catalogue is very satisfac- 
tory, although the width of the CLF for central galaxies is 
significantly broader in the SAM than in the SDSS. To al- 
low for a more meaningful comparison, the lower panels of 
Fig. 3 plot the CLFs obtained from the SAM-GC. The first 
thing to notice is that the width of the CLFs of the central 
galaxies are now in much better agreement with the SDSS; 
apparently, the group finder artificially 'narrows' the scat- 
ter in the relation between halo mass and the luminosity of 
the central galaxy. This simply owes to the fact that we use 
the group luminosity to determine the group mass. Other 
than that, the agreement between the CLFs obtained from 
the SAM-GC, and those extracted directly from the SAM is 
very good, indicating that our group finder allows an accu- 



rate recovery of the true $(L|M) (see also tests in Yang et 
al. 2005b). 

Except for the highest mass bin, the CLFs extracted 
from the SDSS and the SAM group catalogues arc in good 
agreement with each other, indicating that the SAM not 
only fits the gala^xy luminosity function, but it even does 
so as function of halo mass. However, at the massive end 
the SAM predicts significantly more relatively faint galax- 
ies in massive haloes then observed. In the highest mass bin 
shown, the SAM overprcdicts the number of faint satellites 
with L = 3 X 10^ L© by a factor ~ 2. Since these are vir- 
tually all red, early-type galaxies (see below), this suggests 
that the SAM overprcdicts the mimber density of faint, red 
galaxies. Indeed, as already shown in COG, the SAM over- 
predicts the luminosity function of red galaxies at the faint 
end. The analysis here suggests that this largely owes to an 
overabundance of satellite galaxies in massive haloes with 
M > 3 X W^^h-^ Mq. 

We emphasize that this discrepancy is not due to the 
fact that we have ignored fiber collisions in the SDSS. Since 
the spectroscopic fibers of the SDSS have a minimum an- 
gular separation of 55", the spectroscopic catalogue suffers 
from an incompleteness on small angular scales. This will 
impact on the multiplicity function of the groups in our cat- 
alogue. However, as shown by Berlind et al. (2006), the ef- 
fect is relatively small, typically reducing the rrmltiplicity of 
groups by ^ 10 percent, which is negligible compared to the 
factor 2 eluded to above. 

4.2 Blue fraction as Function of Luminosity 

The dotted line in the upper-left panel of Fig. 4 shows the 
fraction of blue galaxies, /biuo, in the SDSS as a function of 
luminosity. Here all 184, 425 galaxies in our SDSS sample de- 
fined in Section 2.1 are used. As is well known, the fraction of 
blue galaxies decreases drastically with increasing luminos- 
ity, dropping from ~ 95 percent at '^'^Mr — 51og/t = —16 to 
< 5 percent for galaxies with " '^Mr — 51og/i < —22.5. The 
dashed line in the lower-left panel shows the blue fraction for 
the 428, 013 model galaxies in the SAM-RS. Although this 
blue fraction also reveals an overall decrease with increasing 
luminosity, there are two marked differences with respect to 
the SDSS. First of all, at the bright end there is a sudden 
upturn in /blue; galaxies with °'^Mr — 5 log ft ~ —23 have a 
blue fraction of ~ 26±6 percent, compared to zero percent in 
the SDSS (note though that there are only 8 SDSS galaxies 
in this luminosity bin). As already mentioned in Section 3, 
these bright blue galaxies in the SAM are ULIRGs for which 
the dust modeling is inadequate (cf. lower panels of Fig. 1). 
The second discrepancy between SAM-RS and SDSS is more 
important; at the faint end the blue fraction in the SAM- 
RS never exceeds 62 percent, and is therewith much lower 
than the blue fraction of faint SDSS galaxies. Consistent 
with what we inferred above from the CLFs, this indicates 
that the SAM severely overpredicts the fraction of faint, red 
galaxies. 

To investigate whether this mainly concerns central 
galaxies, satellite galaxies, or both, we now resort to the 
group catalogues extracted from the SDSS and SAM. The 
open circles with errorbars, connected by solid lines, indicate 
the blue fractions of galaxies that make it into the group 
catalogue. Comparing these for the SDSS to those obtained 
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Figure 4. The luminosity dependence of blue galaxy fractions. From left to right, the panels show the blue fractions of all galaxies 
(centrals plus satellites), central galaxies, and satellite galaxies. The top panels show the results from the SDSS group catalogue (open 
circles with Poissonian errorbars). The dotted line in the top-left panel shows the results obtained from the full SDSS, including those 
galaxies that were not assigned to a galaxy group by the group finder. The bottom panels show the results from the SAM of C06. Results 
are shown for both the SAM-RS (dashed lines) and the SAM-GC (open circles with Poissonian errorbars). See Section 4.2 for a detailed 
discussion. 



from the full sample (dotted lines), we see that the blue frac- 
tion has become somewhat lower at the faint end. We can 
understand this by looking at the blue fractions of central 
galaxies (upper middle panel) and satellite galaxies (upper 
right panel). This shows that both blue fractions decrease 
with increasing luminosity, but that the luminosity depen- 
dence is more pronounced for the centrals. Since the galax- 
ies that do not make it into the group catalogue are mainly 
isolated, and thus central, galaxies that live in haloes below 
the mass limit of the group catalogue, these are mainly blue. 
This explains why the fraction of faint blue galaxies is lower 
in the group catalogue than in the full redshift catalogue. 

In the SAM, the group selection also causes a drop of 
the blue fraction of faint galaxies, but of a much larger 
amplitude. In fact, the fraction of 'group' galaxies with 
-16 > Mb J - Slog /i > -18 is ~ 20 percent in the SAM, 
compared to ^ 60 percent in the SDSS. The reason for this 
discrepancy is largely due to the satellite galaxies: as shown 
in the lower right-hand panel, the blue fraction of satellite 
galaxies in the SAM is much too low, especially at the faint 
end. In the case of the central galaxies (middle panels), the 
agreement between SAM and SDSS is much better, espe- 
cially for centrals with Mtj — Slog/i ~ —21. However, at the 



faint and bright ends, the SAM significantly overpredicts the 
blue fractions by ~ 15 and ~ 25 percent, respectively. 

The dashed lines in the lower middle and lower right- 
hand panels show the blue fractions of centrals and satellites 
of all galaxies in the SAM-RS (including those that are not 
in the group catalogue). This shows that the group finder 
very accurately recovers the blue fraction of central galaxies, 
but slightly overpredicts that of satellites. This owes to the 
interlopers (group members that do not actually belong to 
the same halo), which tend to be isolated, central galaxies in 
low mass haloes, and which are thus preferentially blue. The 
contamination, however, is sufficiently small that it does not 
significantly affect any of our results. 

In summary, although the SAM matches the overall 
blue fraction of galaxies almost exactly (see Section 3.1), 
when split according to luminosity or according to centrals 
and satellites, there are very significant differences between 
SAM and SDSS. The SAM overpredicts the blue fraction 
of central galaxies at both the bright and the faint end of 
the distribution, and dramatically underpredicts the blue 
fraction of (faint) satellite galaxies. In particular, the SAM 
predicts that virtually all ( >, 85%) satellite galaxies are red, 
whereas the SDSS indicates that the fraction of red satel- 
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Figure 5. Blue galaxy fractions as a function of absolute magnitude in the "'^r-band. Results are shown for all galaxies (left-hand 
panels), central galaxies (middle panels) and satellite galaxies (right-hand panels), and for 6 different mass bins as indicated (the values 
in square brackets indicate the range of log(M/h~^ Mq)). Results are only shown for mass-luminosity bins that contain at least 50 
galaxies in total, and for clarity (Poissonian) errorbars are only shown for one mass bin. From top to bottom, results are shown for the 
SDSS group catalogue, the SAM redshift survey (SAM-RS), and the SAM group catalogue (SAM-GC). Note the poor agreement between 
SDSS and SAM-GC, indicating that the SAM is not correctly treating the physics responsible for determining whether a galaxy is red 
or blue. See Section 4.3 for a detailed discussion. 



lites decreases from ~ 90 percent at "'^M^ = —22 -I- 51og/i 
to ~ 40 percent at "'^Afr ~ — 17 + 51og/i. This suggests 
shortcomings for the star formation truncation mechanisms 
in the SAM: apparently the treatment of strangulation is too 
efficient, while the model for AGN feedback is not efficient 
enough (see Section 5 for a more detailed discussion). 



4.3 Blue Fraction as Function of Halo Mass 

We now turn to the mass dependence of /biuc- We split the 
SDSS and SAM group catalogues in six logarithmic mass 
bins and determine how the blue fractions in each of these 
bins depend on luminosity. For each bin in mass and lu- 
minosity, the blue fraction is defined as the total number 
of blue galaxies in that bin, divided by the total number of 
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galcixies in that bin (i.e., wc do not average the blue fraction 
over individual groups or haloes). 

The results for the SDSS group catalogue are shown in 
the upper panels of Fig. 5 and are listed in tabular form 
in Appendix A. The upper left-hand panel shows the re- 
sult for all gala^xies (centrals plus satellites). In each mass 
bin, the blue fraction decreases with increasing luminos- 
ity, but only very mildly. In fact, over the magnitude range 
— 19 >j "'^Mr — 5 log ft >j — 21.5 the luminosity dependence 
is remarkably weak, for all six mass bins. At fixed lumi- 
nosity, however, there is a clear mass dependence, with the 
blue fraction decreasing with increasing halo mass. Over 
the range lO^^/i^^M© < M < lO^'' ft" ^Mq, the blue fraction 
changes by ~ 30 percent, at all luminosities. This indicates 
that the colour of a galaxy is more strongly determined by 
the mass of the halo in which it resides than by its own 
luminosity (cf. Yang et al. 2005b; paper I). Consequently, 
the strong luminosity dependence of /biuo soon in Fig. 4 is 
mainly a reflection of the fact that more luminous galaxies 
typically reside in more massive haloes. Note that some ear- 
lier work (e.g. Balogh et al. 2004, Tanaka et al. 2004) has 
found no correlation between galaxy properties and halo ve- 
locity dispersion. However, as shown in paper I this is most 
likely due to the smaller sample size and the fact that using 
velocity dispersion as a mass estimator tends to smear out 
the mass dependence. 

The upper middle panel shows the blue fractions of cen- 
tral galaxies. Note that at a given halo mass, there is only 
a relatively small dynamic range of luminosities over which 
we can measure fh\ue{L). Nevertheless, at a given halo mass 
there is a clear indication that the blue fraction decreases 
with increasing luminosity. At the same time, at a given 
luminosity the blue fraction also decreases with increasing 
halo mass. Finally, the upper right-hand panel shows the re- 
sults for the satellite galaxies in the SDSS group catalogue, 
which look similar to those for the full sample of galaxies 
shown in the upper left-hand panel. 

The middle and lower rows of panels in Fig. 5 show the 
results obtained from the SAM-RS and SAM-GC, respec- 
tively. In the former no group finder has been applied, and 
the mass binning is according to the true halo masses of the 
galaxies. A comparison with the results obtained from the 
SAM group catalogue is therefore indicative of the accu- 
racy with which our group finder allows a recovery of the 
true underlying /biue(i, M). A number of difi^erences are 
clearly apparent, which mainly owe to the impact of inter- 
lopers and, more importantly, to errors in the assigned group 
masses. Since the true /biue(i, M) has extremely steep gra- 
dients in both L and M, even small errors in any of these 
quantities can have a significant impact on the blue frac- 
tion obtained from the group catalogue. This mainly causes 
errors in the absolute values of /biue(i, M). However, the 
relative relations of fhiue{L\M) and fhiue{M\L) are well re- 
covered. Note that the blue fraction of central galaxies with 
"■^Mr — 5 log ft — —19.5 is less than 10 percent for all mass 
bins shown. Yet, as can be seen from Fig.4, the blue frac- 
tion of central galaxies in the group catalogue with that 
luminosity is ~ 80 percent. This indicates that virtually all 
blue centrals with '^'^ Mr — 5 log ft < —19.5 reside in haloes 
with M < lO^^ft"^ Mq. 

Comparing the fbiue{L, M) obtained from the SAM-GC 
with those obtained from the SDSS, one notices immedi- 



ately that they have very little in common. Probably the 
most dramatic difference between the SDSS and the SAM- 
GC concerns the blue fraction of satellite galaxies (shown in 
the panels on the right hand side), which is much too low in 
the SAM-GC, especially for faint galaxies and for low mass 
haloes. A comparison with the SAM-RS results shows that 
this discrepancy can not be attributed to artefacts of the 
group finder. In addition, the SAM-GC predicts that the 
blue fraction of central galaxies increases with increasing 
luminosity, opposite to what is seen in the SDSS (middle 
panels). This effect is most severe for haloes with masses 
M < lO^^ft"^ Mq. These two problems conspire to pro- 
duce a blue fraction of the full galaxy population (centrals 
and satellites) in the SAM-GC which is very different from 
the SDSS (left hand panels), both quantitatively (mainly 
because of the too low number of blue satellites) and qual- 
ititatively (mainly because of the reversed relation between 
luminosity and the blue fraction of central galaxies). 

We are thus led to conclude that although the SAM 
reproduces the overall blue fraction (when integrated over 
all galaxies), when it comes to fbiue{L, M), there axe dra- 
matic differences between model and data. Note that the 
SAM of C06 fits the overall luminosity function and even 
yields CLFs, i.e., luminosity functions as function of halo 
mass, whose shapes are in remarkably good agreement with 
the SDSS data. This suggests that one of the main problems 
for the SAM is the treatment of the physics responsible for 
determining whether a galaxy is red or blue. This includes 
the star formation truncation mechanisms, such as strangu- 
lation and AGN feedback, but also the treatment of dust 
extinction. In the next section we present a more detailed 
discussion of the possible implications. 



5 DISCUSSION 

The above analysis of the fractions of blue galaxies as func- 
tion of halo mass and luminosity has revealed several prob- 
lems for the SAM of COG. In the following we discuss the 
possible implications for galaxy formation models. 



5.1 Tidal Stripping 

The first problem concerns the abundance of satellite galax- 
ies. As shown in Section 4.1, the SAM overpredicts the num- 
ber of faint satellite galaxies in massive haloes by up to a fac- 
tor ~ 2. Most likely, this indicates that the luminosity evolu- 
tion of the satellites is not properly accounted for. Satellite 
galaxies can become fainter due to star formation trunca- 
tion followed by passive evolution, or due to tidal stripping 
of their stellar mass. Since most of the satellite galaxies in 
the SAM are already red, they can not become much fainter 
than they already are for their given stellar mass. In other 
words, there is little to gain from adding physical processes 
that may speed up the star formation truncation, such as 
ram pressure stripping. In fact, this will only increase the 
fraction of red satellites, which is already much too large. 
The most plausible explanation for the overabundance of 
satellite galaxies is the neglect of tidal stripping. It is well 
known that the tides in massive haloes can easily strip satel- 
lite galaxies and their dark matter haloes of large fractions of 
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their mass. This is supported by the detection of intraclus- 
ter light (e.g., Bernstein et al. 1995; Gonzalez et al. 2000) 
due to a diffuse population of intergalactic stars. Most likely, 
these stars, which contribute around 10 percent of the total 
cluster light (Zibctti ct al. 2005), have been tidally stripped 
from satellite gala^xies. 

Formalisms to describe tidal stripping in the presence of 
dynamical friction in a background potential have been de- 
veloped by, among others, Taylor & Babul (2001, 2004) and 
Zontncr & Bullock (2003). As shown in Benson ct al. (2002), 
including such a formalism in semi-analytical models signif- 
icantly reduces the abundance of satellite galaxies. Further- 
more, since tidal stripping predominantly affects satellites 
with small pericentric radii, which arc typically redder (blue 
galaxies have only recently been accreted by the halo, and 
have not yet experienced much dynamical friction), tidal 
stripping may also reduce the red fraction of satellites and 
thus help towards solving the problem with the red fractions 
of satellites being too large. 

5.2 Strangulation 

The second problem for the SAM of C06 is that the fraction 

of blue satellite galaxies is much too low, especially for faint 
galaxies and in low mass haloes. This suggests that 'stran- 
gulation', as incorporated in the SAM, is much too efficient. 
In virtually all semi-analytical models, strangulation is in- 
cluded and modeled in the same way: as soon as a galaxy 
becomes a satellite galaxy, its hot gas reservoir is 'stripped' 
(i.e., the hot gas that belonged to the satellite galaxy is 
added to the hot gas reservoir of the central galaxy). Conse- 
quently, after a delay time in which the new satellite galaxy 
consumes (part of) its cold gas, its star formation is trun- 
cated. The results presented here suggest that this formula- 
tion is too crude, as it predicts a blue satellite fraction that is 
much too low^. In particular, strangulation is modeled with- 
out any explicit halo mass dependence, which explains why 
the blue fraction of satellite galaxies in the SAM is virtually 
independent of halo mass. In the SDSS, however, the blue 
satellite fraction decreases with increasing halo mass, sug- 
gesting a clear mass scaling of the strangulation efficiency. 
The physical mechanisms thought to be responsible for the 
removal of the satellite's hot gas reservoir are tides and ram 
pressure stripping. The latter requires that the parent halo 
has a sufHciently dense hot corona. Since this requirement 
is more likely to be fulfilled for more massive haloes, which 
in general have a larger fraction of hot gas, one actually ex- 
pects that the strangulation efficiency increases with increas- 
ing halo mass. Indeed, using numerical simulations, Bekki et 
al. (2002) find that strangulation is significantly more effec- 
tive in massive galaxy clusters than in lower mass groups. 
It remains to be seen whether a simple scaling along these 
lines can bring the red and blue fractions of satellite galaxies 
in SAMs in better agreement with the SDSS data presented 
here. 

Interestingly, in a comparison of semi-analytical models 
and hydrodynamical SPH simulations of galaxy formation, 

^ Note that the SAM of C06, as most other SAMS, does not even 
account for ram pressure stripping, which will only shorten the 
truncation time, thus decreasing the blue satellite frsLction. 



Zheng ct al. (2005) have shown that the later predicts that 
haloes with 10^^ Mq < M < lO" Mq have a significantly 
higher fraction of young satellites (similar to the blue satel- 
lites discussed here) than SAMs, in much better agreement 
with the SDSS results presented here (see their Fig. 4). Sim- 
ilarly, Cattaneo et al. (2006b) find that their semi-analytic 
model produces too few blue satellite galaxies compared to 
their SPH simulation. Since these SPH simulations foUow 
the actual hydrodynamical processes leading to strangula- 
tion, this indeed suggests that a more realistic treatment of 
strangulation in the SAMs may solve the problem with the 
colors of satellite galaxies. SPH simulations such as those 
described in Zheng ct al. (2005) and Cattaneo et al. (2006b) 
may prove useful in calibrating such a new and improved 
strangulation model. 



5.3 Dust Extinction and AGN Feedback 

The SAM also has problems with the blue fractions of cen- 
tral galaxies, which are too high, especially at the bright 
and faint ends. In addition, for a given halo mass, the blue 
fraction of central galaxies increases with luminosity, con- 
trary to what is seen in the SDSS. Which aspect(s) of the 
semi-analytical model are responsible for these problems is 
not entirely clear. They can indicate a problem with the 
modeling of dust extinction, a problem with the treatment 
of AGN feedback, or both. The former is almost certainly 
responsible for the overproduction of bright and blue cen- 
trals. As already discussed in C06, this population of galax- 
ies is reminiscent of the ULIRG population, for which the 
oversimplified treatment of the dust extinction is certainly 
inadequate: real starburst are likely to be accompanied by 
additional extinction which would make these galaxies both 
fainter and redder. This would help to suppress the strong 
increase of /blue with increasing luminosity, though it re- 
mains to be seen whether it can result in a blue fraction that 
decreases with increasing luminosity, as observed. Further- 
more, in order to suppress the impact of these starbursting 
model galaxies on the comparison presented here, we already 
removed all galaxies with {g — r) < from the SAM. De- 
spite this, however, the SAM still significantly overpredicts 
the fraction of blue centrals with "'^Mr — 51og/i ~ —23. 
Although improper dust modeling is likely to be the cause 
for this discrepancy, it remains to be seen whether this can 
also explain the overprediction (by ~ 15 percent) of the blue 
fraction of faint centrals. 

Alternatively, the blue fractions of central galaxies may 
be modified by changing the AGN feedback description. As 
discussed in the introduction, it is still largely unknown how 
AGN impact on their surroundings, and thus feed back on 
the process of galaxy formation. It should therefore not come 
as a surprise if the parameterization of C06 is not entirely 
correct or complete. The purpose of this paper, however, 
is not to find an improved formulation of AGN feedback. 
Rather, we have presented data, in the form of blue galaxy 
fractions as function of both halo mass and luminosity, which 
we believe to be useful in discriminating between different 
AGN feedback models. Future SAMs can test and/or cali- 
brate their particular parameterizations against these data. 
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6 CONCLUSIONS 

It has become clear that a successful reproduction of the 
galaxy luminosity function, and of the bimodality in the 
galaxy distribution, requires a mechanism that can trun- 
cate the star formation in massive galaxies at relatively late 
epochs. At the same time, the fact that red, passive galaxies 
preferentially reside in overdense regions, such as clusters 
and groups of galaxies, suggests that star formation trun- 
cation acts preferentially in massive haloes. Current mod- 
els consider two such truncation mechanisms: strangulation, 
which acts on satellite galaxies, and AGN feedback, which 
predominantly affects central galaxies. 

Typically, galaxy formation models axe tuned to repro- 
duce the global properties of the galaxy distribution, such as 
the luminosity function and the total fraction of blue and 
red galaxies. However, oven when two models predict exactly 
the same global statistics, their statistics as function of halo 
mass may be very different. The latter is clearly more con- 
straining for the model, and furthermore, holds important 
information regarding the mass-dependence of the various 
physical mechanisms associated with galaxy formation. In 
particular, since star formation truncation causes a galaxy 
to become red, the relative fractions of red and blue galaxies 
as a function of halo mass holds important clues regarding 
the halo mass dependence of the efficiencies of AGN feed- 
back and strangulation. 

To provide a test-bed for models of galaxy formation, 
we have used a galaxy group catalogue extracted from the 
SDSS for which we have computed the fraction of blue and 
red galaxies as function of both galaxy luminosity and group 
(halo) mass. To illustrate the potential constraining power 
of this data we have compared these fractions to those in 
the semi-analytical model for galaxy formation of Croton et 
al. (2006), which includes both 'radio mode' AGN feedback 
and strangulation. To allow for a fair comparison between 
the SDSS group catalogue and the SAM, not influenced by 
potential inaccuracies associated with the group finder (i.e., 
interlopers, incompleteness, errors in assigned group mass), 
wo have applied the same group finder over a mock rcdshift 
survey constructed from the SAM. The fbine{L, M) obtained 
from this SAM group catalogue is in fair agreement with the 
true /biue(i, M) obtained directly from the SAM, indicating 
that our group finder allows a reliable recovery of the blue 
fraction as a function of both galaxy luminosity and halo 
mass. Although interlopers and errors in the group masses 
may cause some errors in the absolute values of /biue(^, M), 
the relative scalings of f\,i-ae{L\M) and /biue(Af|L) are well 
recovered. 

Although the SAM fits the overall luminosity func- 
tion, reproduces the overall color distribution of the SDSS 
galaxies, and even predicts a conditional luminosity func- 
tion whose shape is in excellent agreement with the data, 
its prediction of /biue(i, is in poor agreement with the 
SDSS data. In particular, we have identified four problems: 

(i) In massive haloes the abundance of faint satellite 
galaxies is too high by up to a factor ~ 2. 

(ii) The fraction of blue satellite galaxies is much too low, 
especially for faint galaxies and in low mass haloes. 

(iii) The fraction of blue central galaxies is too high, es- 
pecially at the bright and faint ends. 

(iv) For a given halo mass the blue fraction of central 



galaxies increases with luminosity, contrary to what is seen 
in the SDSS. 

The first of these problems is likely to owe to the fact that 
the SAM does not model the tidal stripping of the stellar 
mass of the satellite galaxies as they orbit the parent halo. 
As shown in Benson ct al. (2002), inclusion of this effect 
significantly reduces the abundance of satellite galaxies at 
a fixed luminosity. The second problem is most likely due 
to an oversimplification of the treatment of strangulation. 
In the SAM, strangulation occurs instantaneously, indepen- 
dent of halo mass. However, based on the SDSS data, we 
have argued that the strangulation efficiency has to scale 
with halo mass, such that more massive haloes strangulate 
their satellites on a shorter time scale. The physical motiva- 
tion for such a scaling is that the ram pressure stripping of 
the hot gas reservoir of newly accreted satellites requires the 
parent halo to have a sufficiently dense corona of hot gas. 
Since the fraction of hot gas is typically an increasing func- 
tion of halo mass, this may introduce a mass dependence 
in the strangulation efficiency as required. Indeed, hydrody- 
namical SPH simulations, which automatically take this into 
account, seem to predict blue satellite fractions that are sig- 
nificantly higher than in the semi-analytical models (Zheng 
et al. 2005, Cattanco ct al. 2006b). Finally, the third and 
fourth problem listed above, both of which concern central 
galaxies, arc likely to reflect shortcomings of the modeling 
of dust extinction and/or AGN feedback. 

In summary, galaxy formation models are often tested 
and calibrated against global properties of the observed 
galaxy distribution. Recently, with the inclusion of AGN 
feedback, numerous studies have claimed success in repro- 
ducing these global statistics, even though very different for- 
mulations for the various physical processes have been used. 
In order to discriminate between these models more spe- 
cific data is required. In this paper we have presented the 
fractions of blue and red galaxies as function of luminosity, 
halo mass, and separately for central and satellite galaxies. 
Clearly this data is far more constraining, and thus more 
challenging, than the global fractions of red and blue galax- 
ies, or the overall luminosity function. We have shown that 
indeed these data provide valuable new insights regarding 
the physics of galaxy formation, and we hope that they will 
provide a useful test-bed for future models of galaxy formar 
tion. 
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APPENDIX A: BLUE GALAXY FRACTIONS IN 
THE SDSS 

The following three tables list the fraction of blue galax- 
ies for all galaxies (Table Al), for satellite galaxies (Tar 
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Table Al. Blue fractions of all galaxies 
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Table A2. Blue fractions of satellite galaxies 
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ble A2), and for central gala^xies (Table A3), as obtained 
from our SDSS group catalogue. Columns correspond to dif- 
ferent bins in log(M), with M in /i~^Mq, as indicated at the 
top in square brackets. Rows correspond to different magni- 
tude bins ("'^Mr — 5 log ft), as indicated at the left in square 
brackets. Each entry lists the blue fraction plus, in brackets, 
the total number of galaxies (all, satellite or central) in that 
bin. As in Fig. 5, only entries with at least 50 galaxies are 
indicated. 
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Table A3. Blue fractions of central galaxies 



[12,12.5] [12.5,13] [13,13.5] [13.5,14] [14,14.5] [14.5,15] 
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